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Abstract. Recently the block-diagonalization of Dirac-operators was inves- 
tigated from a mathematical point of view in the one-particle case 1131 . We 
extend this result to the TV-particle case. This leads to a perturbative realiza- 
tion of the Furry picture in the TV-particle two-spinor space. 



1. Introduction 



The idea of block-diagonalizing the Dirac operator, i.e., decoupling electronic and 
positronic states in such a way that the upper components of a 4-spinor correspond 
to electronic and the lower components to positronic states, goes back to Foldy and 
Wouthuysen They succeeded in decoupling the free Dirac operator and also 
addressed the case with interaction in the non-relativistic limit. Unfortunately, their 
expansion does not converge (see Thaller chapter 6] and references therein), 
j^. . The reason is that the correct parameter for the expansion, in order to have 

convergence of the spectrum, is not the inverse velocity of light but the coupling 
constant of the external potential. A perturbative and iterative method to accom- 
plish this is due to Douglas and Kroll 3 and was corrected by Jansen and Hcfi 
^ '. 0. The method is very attractive for numerical calculations because the resulting 

Oh' Hamiltonians operate on two-spinors and has been successfully used for calculations 

^ (— i in relativistic quantum chemistry in the last twenty years (see El EU 

and references therein). 

From a mathematical point of view, the one-particle case was investigated re- 
cently by Siedentop and Stockmeyer ^3] (see also ^2]). They proved, under suitable 
conditions on the potential, that there exists a family t/ 7 of unitary operators, an- 
alytic in the coupling constant 7, that exactly decouples electronic and positronic 
states. Moreover, it was shown that the block-diagonalized Dirac operator devel- 
oped in power series in 7 coincides, at least formally in the first orders, with the 
operators resulting from the method of Douglas, Kroll and Hefi. Moreover, they 
proved that the spectra of the truncated expansions converge to the spectra of the 
original operator. 

The aim of this work is to extend their results to the iV-particle Coulomb-Dirac 
Hamiltonian. We consider the TV-particle Coulomb-Dirac operator in the Furry 
picture, i.e., the operator restricted to the positive spectral subspaces of each one 
particle operator. Using simple generalizations of the methods used in |14| we prove 
norm resolvent convergence of the operators resulting from the power expansion of 
the projected Hamiltonian in powers of the coupling constant. For convenience for 
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the reader, we give in section [5] the main definitions and results of |14j . as far as we 
need them. 

2. Definition of the problem and notation 
The N-particle Hilbert space is denoted by 

#W =3®... ®j5, 

where Sj := L 2 (R 3 ,C 4 ). An extension of a closable operator A on f) with domain 
V to $j( NS> acting on the j-th component is written as 

Aj := 1 ® ••• ® <g> ••• <g) 1, 

j-th place 

and the n-fold tensor product as 

A:= A® ■■ ■ ® A. 

If A is essentially self-adjoint on V, then all of the above operators are essentially 
self-adjoint on ® J _ 1 I> (see Jjjj, chapter VIII. 10). For technical reasons, we will 

need the operator 2)$ :— \Do\j, where Dq is the free Dirac operator. 

The Coulomb-Dirac operator is given formally by 

N 

H^--=Y,[D^ + l E W*> (!) 

j=l l<i<j<N 

where, using units h = c = 1, Z is the atomic number of the nucleus, 7 = Za 2 is the 
coupling constant of the one-particle potential, and a is the fine structure constant. 
D 1 = Do + -fV where Dq is the free Dirac operator, V = — 1/| • | is the Coulomb 
interaction with the nucleus, and Wij is the interaction between the particle i and 
j, which is defined by 

{W lj f){x 1 ,...,x N ) := ; — - — f (2) 

\ x i x j\ 

for / € i/ 1 (M 3 ) 4 <gi ■ ■ ■ ® H 1 (R 3 ) i . We set W 1 := § Ei< J<j <jv W H- lt is wel1 
known that the operator in Q (without the electron-electron interaction W 7 ) has 
the whole real axis as its spectrum. 

We consider instead of the Coulomb-Dirac operator in the Furry picture, that 
is we restrict Q onto the positive spectral subspaces of each one-particle operator 
D 1 . The Hilbert space is given by 

where Pj. :— X[o,oo) {D-i)- The TV-particle projection is given by 

VI :=Pl®---®Pl. (3) 

We are interested in VIHq^VI in particular, in its realization in the N particle 
two-spinor space. We recall the definition and properties of the unitary transforma- 
tion Uj given in ^1] (see equation 1)22(1 below) and consider its TV-particle version 

U 7 := Z7 7 (g) • • • (g) U T (4) 
This operator has the property U-fPllA^ 1 = V\. Moreover, the TV-particle Foldy- 
Wouthuysen transformation Ufw fulfills UfwP+ = P+ ® • • • ® P+Ufw where fi+ :— 
(1 + f3)/2 is the projection onto the upper two-spinor. Therefore, the formal Hamil- 
tonian in the Furry picture realized in the TV-particle two-spinor space is given by 

Hg a& = u F ^u 1 viH^2'Plu; 1 u^. (5) 
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Our main result, Theorem [5] below, is that the spectra of the operators resulting 
from the expansion in 7 of J^J) converge to the spectrum of the Furry-operator given 
formally by V\H^V\. 

3. Main result 

The one-particle Dirac operator with Coulomb potential is given by 

zV=a-iv+/3-7rq 

acting in the Hilbert space L 2 (R) 4 . For |-y| < \/3/2 the operator D 7 is essentially 
self-adjoint on C£°(R 3 \ {0}) 4 and self-adjoint on H 1 (M 3 ) 4 . Furthermore, for 7 < 1 
the operator has a distinguished self-adjoint extension characterized by D(D 1 ) C 
if 1 / 2 (see Thaller [T31 Theorem 4.4.]). Moreover, for I7I < 1 we have that |D 7 | > 
v/l - V>0. 

We define the operator T 7 : — > f)^ V ' ) without electron-electron interaction 

as 

f, = rlJ2[D^vl, (7) 

i=l 

with domain S3(T 7 ) = VlH^R 3 , C 4 )® N . For | 7 | < y/3/2 the operator f 7 is 
essentially self-adjoint (0 Theorem VIII.33]). We denote the closure of T 7 by T 7 
and its form domain by 0(T 7 ). 

Define the quadratic form q(f,g) := (f,T^g) + (f,W y g) for f,g E 0(T T ). We 
have the following: 

Theorem 1. Let < 7 < V3/2. There exists a unique self-adjoint operator 

fT 7 + : £^(7) -> fl+°(7) £(#7,+) = ^( T 7) = Pj-H" 1/2 (M 3Ar ,C 4N ) w/iere 
^(iJ^R 3 ) 4 )®^ is a /orm core /or iJ 7i+7 such that for all f,g G 0(# 7) +) 

5) = {f,H 7 ,+g). 

This self-adjoint extension is the Friedrich extension of the symmetric operator 
T 7 + V\W^P\ defined on V\ ®f =1 i? 1 (R 3 ) 4 . 

Proof. Step 1: Well-definedness of q. We have by LemmaE! 

(/,|A>kf) < j-UAD^if) < i(/,r 7 /). (8) 

This inequality together with Lemma ^ implies the well-definedness of q. 
Step 2: Definition of H 1:+ . We mimick the proof of the KLMN-theorem, using 
the notation of of Reed-Simon [THl Theorem X.17]. Pick / e 0(T 7 ). We start by 
proving that (/, W 7 /) < c(/, T 7 /) for some c > 0. Since for / G 55+ , we have 

= 1 E (^^Vj^l^l^oir 172 ^^!^!- 172 !^!^ 1 / 2 /) 

l<i<j<N (9) 
^ 7?r W-l/2f ID I T-l/2fl / 7£iV^iV--l) 2 2 

- 2Z ( T 7 /jPoli-Ty 7 /) < ^ 11/11 =:c||/|| , 

l<i<j<N 7 

where we used Lemma ^ and equation (j3J). Therefore, 

(/, T,f) + (/, W,f) < (1 + c)(/, T 7 /) < (1 + c)(/, (T 7 + W 7 )/). (10) 

The latter shows that the norms || • ||+i,t„ and || • are equivalent. Thus q is 

a semi-bounded, closed quadratic form on 0(T 7 ), which therefore defines a self- 
adjoint operator -ff 7 ,+ with form-domain 0(iJ 7i+ ) = 0(T 7 ). 
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Step 3: Determination of the form domain. We have the following chain of in- 
equalities, where we used Lemma El in the second and Hardy's inequality in the 
third step: 



N N 

i=l i=l 
N 

5^(1 + 2 7 )|A,|i < (1 + 2j)Ny/-A 3N + 1 (11) 
t=i 

Inequalities (|ll|l imply that if 1 / 2 (M 37V , C 4 ) is complete with respect to the qua- 
dratic form b(f,f) := (ySiT^kf, ySTSW) for / G H 1 /*^, c ^ } . 
Since ^ commutes with \D^, we have 7>7_ # ^ 2 (R 3N , C 4 ™) C tf 1 / 2 ^, C 4 ™). 

Moreover, ^i/ 1/2 (]R 3Ar , C 4 ™) is a closed subspace of H^ 2 (R 3N , C 4 ") with respect 
to the norm generated by b. Hence 'P7iJ 1 / 2 (lR 3A ', C 4 ™) is complete with respect 
to the restriction of b. Since T 7 is essentialy self-adjoint the self-adjoint operator 
associated to the restriction of b is T 7 . Thus, 0(T 7 ) = VlH^ 2 {M. 3N ', C 4 ™). 

Since obviously 3) C G(^+) = 7^H 1/2 (R 3Ar , C 4 "), we have that 
is the Friedrich extension of the symmetric operator T 7 + P7W 7 'P7 defined on 

We now turn to the operator 1l 1 := ® ■ ■ ■ ® U 1 which is a unitary map- 
ping U-f : .$3+^(7) — > Sj^\o). We define the operator -ff 7 ,+ : Sj^q — > Sj±£ 
as := U y H^ t +U~ 1 , where : ~ #+^(0). Analogously, interpreting the 

Foldy-Wouthuysen transformation in a natural way as a mapping Wfw : $4^0 ~~ * 
L 2 (R 3 , C 2 )® N = L 2 (R 3A \C 2 ™), we define : L 2 (R 3N , C 2 ™ ) -» L 2 (R 3A \C 2 ™) 

as the operator i? 7i + s = Wfw#7,+^fw- I 11 wa y> -^ts+ S can ^ e seen as the block- 
diagonalization of ff 7 ,+- By Theorem^ i?-y,+ and if^ 1 !]? are self-adjoint operators 
with form domain W 7 Q(T 7 ) and WfwW 7 £}(T 7 ) respectively. Actually £/ 7 0(T 7 ) = 
V%H l / 2 (R m , C 4 ™), since 5> 1/2 W 7 ^~ 1/2 is a bounded operator (see proof of Lemma 
EJ. Since W FW commutes with O , we get Wp W W 7 £l(T 7 ) = J ff 1 / 2 (R 3Ar , C 2 ™). 

We denote by hJl + the formal Taylor expansion of H^ 13 ^ up to the power 7^ 

inclusive, acting on C := UFw^-y'P+ (l)H 1 (R 3 , C 4 )® N which is a form core for i? 7 .+ s . 

We set i2* = - fr* + . 

The main result of this letter is: 

Theorem 2. There exists a j c > such that for < j < j c the operators + 
admit a distinguished self-adjoint extension h£ + for k big enough with the property 

that D(/i 7)+ ) C Q(Hf; + g ) = H 1 ' 2 ^ ,C 2N )- Moreover + -» H^f as k -» 00 
in i/ie sense of norm resolvent convergence. 

Proof. According to Kato [HI Theorem VI. 3. 11 and Corollary VI. 3. 12] it is enough 
to prove there exist a sequence afc with ctfc — ► as k — > 00, such that for any f E C 

(/,^/)<a fc (/,<T/)- (12) 

This is equivalent to 

(/, H-'fu-'U^RpFwU^H^ff) -> (13) 
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for / e V+i^H^R^C 4 )®" . Then 

< \\%- 1/2 r^o 1/2 \\ H'%% ll2 f H /2 H-}ff\\\ 

the last term goes to zero due to lemmas El El and 0] below. We used in the last 
step that Ufw commutes with \Dq\ and therefore Ufw with @q and any function 
of it. □ 

Remark 1. The unitary transform U y is not unique. If we take the choice given 
in |14) 11' I'D in the appendix, we know that 7 C > 0, 3775 which corresponds to the 
critical atomic number Z = 52. 

4. Auxiliary lemmas 
The following bound is Kato's inequality. 
Lemma 1. Let f G H^ 2 {R 3 N, C 4 ™). Then 

(/,W«/)<!(/.IW)/°r<^j (15) 

Proof. Pick / € ®f =1 i^QR 3 ) 4 arbitrarily. 



(f,Wijf) = I dxi - ■■dx N f(x 1 , ...,Xi,.. .,x N )- — - — r/(xi, ...,Xi,... ,x N ) 

\Xi Xj\ 



dxi ■ ■■dx N f{xi,. ..,Xi+Xj,.. .,x N )-. — :f(xi, ...,Xi+Xj,.. .,x N ) < 

\ x i\ 

^ J dxi ■ ■■dx N f(x 1 , ...,Xi + Xj,.. . ,x N )\D \if(xi, . . . ,x t + x j: . . . ,x N ) < 

|(/,po|i/)> (16) 

where we used Kato's inequality. This inequality extends by continuity to / € 
H^ 2 (R 3 N,C 4N ). □ 

Lemma 2. For \j\ < "f c the operator S> 1 ^ 2 H^ / la ^£) is bounded and real analytic 
around zero. In particular \\% 1/2 R*% 1/2 1| -> 0. 

Proof. First note that \Dq\J 2 Sl^ 1 ^ 2 is bounded. To prove this, we take an arbitrary 
/ € <g)f =1 H 1 ^ 3 ) 4 and note that 

IH^oL 1/2 /l| 2 = (/.Poll/) < (/,Po|i/) + ••• + (/, \D \ N f) = H /2 f\\ 2 , (17) 

the claim follows by density of H^R 3 ) 4 in £ (7V) . 

According to LemmaEl the operator iDol -1 ^ 2 ^!^^ Po|~V 2 is bounded and 

analytic in 7 for |7| < j c , and so is [(\D \^ 1/2 U 1 D 1 U*\D \^ 1 ^ 2 )} m for m=l,...,N. 
In order to prove the claim it suffices to show the analyticity of each summand of 

One-particle terms: The analyticity follows immediately from 

% ll2 \U 1 PlD 1 PlU- x \ m % 112 = V\% ll2 W 2 \ m 

x [iDol-^U-^U^Dor^UDoll/ 2 ^- 1 ^. (18) 
Interaction terms: The operator 

po|- 1/2 W 7 Po|m /2 = U 7 <g> • • • ® |A)r 1/2 E/ 7 Po| 1/2 ■ • • ® 
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is analytic by LemmaEl and l^olm^Wm.il-Dolm 1 ^ 2 is bounded by Lemma^ Thus, 
writing 

P^o" 1/ Vo|^ /2 |i5o|^ / XPo|^ 2 |I>o|^ /2 T^,i|^oL 1/a 

\d \u 2 u^\d \-^\d \u 2 ^ 1/2 K, (19) 

we have shown the analyticity of the interaction term. 

Since Wfw commutes with @ the operator S^q 1 ^ 2 H^^&q 1 ^ 2 is also analytic 

and therefore has a convergent Taylor expansion with \\@ 1 ' 2 Bt@Q 1 ^ 2 \\ — > as 
k — > oo. □ 

Lemma 3. For \^\ < j c the operator ^q^ 2 U 7 S> is bounded on 

JV 

Proof. We have that S^ 72 < |P>o|* 72 . Note further that by LemmaODthe oper- 

i=l 

ator |£)o| 1/2 ^ 7 |D r 1/2 is bounded, thus 

X^ - 1/2 n = K /2 (]T i^oi- 72 ) -1 E i^if^o" 



N N 

*o /2 Uy®o 1/2n - na 1 / J ^in.i 1 /^-iV l in 1 -// 1 2 

iV JV 



= ||^ 1/2 (E l^olP)" 1 Ed^ol-^l-Doir^l-Dol- 72 ^" 172 )!! 
»=i i=i 

JV 

i=l 

JV 

^^llpoiyXlDo^lllll^oll 72 ^- 17 -' 1 

i=l 



(20) 

is finite. □ 
Lemma 4. For / G #^(7), andO < 7 < ^3/2 we kw ifte estimate \\^ /2 H^ 1 / 2 f\\ < 

i/V^II/ll- 

Proo/. Pick / G P+M (i^QR 3 ) 4 )®^. Then using LemmaEJand W 7 > we get 



(f,@of) < ±-{f,Y,( P l^Pl)jf) < ^(f,H^, + f). (21) 

□ 



' 3 ' 



5. The one- particle case 

Let us define the following operator 

U y := (P°jRj! + P°P2)(1 - (P° - Pi)) 2 )- 1 / 2 . (22) 

Some important properties Q14I Theorem 1, Theorem 2, Lemma 9]) of Z7 7 are listed 
in the following lemma. We set 7 C := 0.3775. 

Lemma 5. (1) Uj is analytic in^ and unitary for \<y\ < 0.6841 and fulfills the 
relation 

U^Pg =PlUj. 

(2) The operator \D \ 1/2 U^\D \- 1/2 is bounded and analytic in 7 for 7 < j c . 
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(3) The operator \Dq\ 1 / 2 U 1 D 1 U*\Dq\ x I 2 is bounded and analytic in 7 for 
7 < 7c- 

The following inequality is used in this paper. For < 7 < set C 7 := 
I ( V 4 (7) 2 + 9 - 4 7 ) and d. ( := ±(1 + C 2 - yj(l - C 2 ) 2 + ^ 2 C 2 ). We have (see 
Morozov [H] and also Brummelhuis et al 2 ): 

Lemma 6. For < 7 < the operator inequality 

|£> 7 | 2 > 4|A)| 2 (23) 

holds. 
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